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Contrary to the longtime and widely conceived belief, we proved that the specific heat coefficient y -also 
called Sommerfeld coefficient - of the interacting Fermion system is not renormalized by the interaction as far 
as the system remains a Fermi liquid state. 
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Introduction - Fermi liquid theory|lH3[] is the most funda- 
mental conceptual building block of the modern quantum the- 
ory of the interacting fermion systems such as metals, semi- 
conductors, superconductors, liquid 3 He, neutron stars, etc. In 
a nutshell, it suggests that an interacting fermion system can 
be one-to-one mapped to a non-interacting fermion system for 
the low energy excitations. In the process of this adiabatic 
mapping, the essential effect of the interaction is to renormal- 
ize the original bare fermions into a renormalized fermionic 
"quasiparticles". While the charge and spin quantum num- 
bers - when they exist - of the fermions are protected by the 
gauge invarience[4] and relativity, the mass of the fermion in 
the condensed matter is an effective mass from the beginning 
and can be renormalized to be a different value from the orig- 
inal effective mass mo of the non-interacting limit. There- 
fore, the renormalized effective mass m* of the quasiparticle 
(q.p.) is the single most important quantity which determines 
the low energy properties of the interacting fermion systems. 
Hence, the reliable measurement of this quantity by experi- 
ments should be of principal importance to study the nature 
and strength of the interaction of the fermionic system. 

There are several different probes to measure the effec- 
tive mass: specific heat (SH) coefficient, de Haas-van Alphen 
(dHvA) effect, angle resolved photo-emission spectroscopy 
(ARPES), optical spectroscopy, etc. Although some interpre- 
tations might be necessary to extract the value of m* from the 
above listed measurements, theoretically all these measure- 
ments should give consistently the same information about the 
effective mass m*. For example, the ARPES measures the q.p. 
energy dispersion E(k) vs momentum k and directly shows us, 
without interpretation, how heavy or light the q.p.s move. The 
dHvA effect similarly depends on the q.p. dispersion E(k), 
so that its measurement also provides a direct information of 
the renormalized mass. However, since the construction of 
the Landau Fermi liquid phenomenology! 1] and its theoreti- 
cal justifications]^ |j, |5[ |6J], the most commonly used probe 
for the effective mass of the q.p.s in the Fermi liquid systems 
is the measurement of the SH coefficient. In particular, Lut- 
tinger had shown in his seminal paper ||3|] in 1960 that the SH 
coefficient y (= lim^o C(T)/T) should be enhanced from the 
non-interacting value yo such as y/yo = m* /mo. Since then, 
the measurement of y has been established as the most impor- 
tant tool to measure the effective mass of the fermionic q.p.s. 
in the condensed matter systems. 



In this paper, we showed that there is an error in the proof 
of Luttinger and the SH coefficient in the interacting fermion 
system is not renormalized by the interaction, hence y/yo = 1 
regardless of the strength of the fermion-fermion interaction. 
Our finding should have far reaching consequences in the 
study of various interacting fermion systems such as strongly 
correlated metals, liquid 3 He, neutron stars, etc. In this paper, 
we will be focusing only on the questions of where was wrong 
in the Luttinger's proof and what is the correct answer for the 
SH coefficient y of the interacting fermion systems. 

SH coefficient y and DOS - It is well known that the SH 
coefficient of the non-interacting fermion system yo is given 
by 



KmC(T)/T=y = -N (0), 



(1) 



where No(Q) is the density of states (DOS) of the non- 
interacting fermion system at the chemical potential. Intu- 
itively, the SH coefficient of the interacting fermion system 
y is expected to be given with the above equation by replac- 
ing No(0) by the DOS of the interacting fermion system N(Q) 
such as y = ^-N(0). But this absolutely reasonable intuition 
falls in a serious trouble as follows. The DOS N(0) of the 
interacting fermion system can be calculated if we know the 
exact one-particle Green's function which is formally written 
as G(k,(a) = ^z^pm^ with the exact self-energy E(fe,co). 
And it is almost trivial to show that N(Q) = Nq(0), as shown 
below, which leads us to conclude that y = yo- However, this 
conclusion is in stark contrast to the common knowledge that 
the SH coefficient should be enhanced by interaction such as 
y/yo ~ m* /mo > 1. There are two possible options to resolve 

2 

this dilemma: (1) y= \-N(0) is not true for the interacting 
system; or (2) the common belief y/yo ~ m* /mo is wrong. 

Let us begin with showing N(Q) — M)(0). It is well known 
that lim7;<D^o£(£,co) = — ^co — ;8 in the Fermi liquid state 
iSll. Then 
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where the wave function renormalization parameter 
A.£ (Zt > 1) is defined. Therefore, N(0) = No(0) is a rigorous 
identity as far as the Fermi liquid state holds even with an 
arbitrary momentum dependent self-energy £(&, co). 

To demonstrate the correctness of the result of Eq.(5), we 
show the numerical results of iV(co) = — ^ ^/mG(fc, CO ) of a 
simple toy model in Fig. 1(a) with varying strength of inter- 
action. In this examplary calculations, we assumed a box 
like DOS for the non-interacting fermion system as iVo(co) = 
1 .0 for —A < CO < A and the effect of interaction is simu- 
lated by the Fermi liquid type self-energy TmE(co) = aco 2 for 
—A < co < A including the corresponding real part /teE(co). 
We chose A = 5. The results are self-explaining, showing 
N(Q) = Nq(Q) for all interaction strength of a. Increasing the 
interaction strength, the width of the q.p. DOS around co = 
becomes progressively narrowed and the spectral weight out- 
side of it is depleted toward the high energy region which is 
not fully displayed here but the total spectral weight of the 
DOS should be conserved. The width of the q.p. DOS around 
CO = is roughly proportional to ~ 1/Z and the value of Z is 
determined by the combination of the values of the interaction 
strength a and the band width scale A. 

SH coefficient J of Interacting Fermi Systems - To find an 
exact theoretic formula to calculate the SH coefficient y of the 
interacting Fermi systems, we start with the same Hamiltonian 
for the interacting fermion system used by Luttinger and Ward 
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(6) 



where 8 r is the energy measured from the chemical potential 
of the non-interacting single particle states with the index r = 
(k,a) for both momentum and spin. c y r , c r are the creation and 
annihilation fermion operators, respectively, and (r*|v|r s ) is 
the general four point fermion interaction matrix. 

In Ref.lQ], Luttinger and Ward wrote down the celebrated 
free energy functional of the interacting fermion system as 



fl(r) = -r£ 



e to " 0+ {ln[e r + E r (co„)-/co„ 



-G,-(co„)E r (co„)} + n 



(7) 



where C0„ = KT(2n + 1) is Matsubara frequency. G,-(co„) and 
E, (co„) are the full Green's function and the full proper self- 
energy, respectively. The functional Q. is defined by LW 
(refers Ref.[2]) as 
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FIG. 1: (Color online) (a) The full DOS AT(co) = Z k ImG{k,(a) 
with a phenomenological Fermi liquid type self-energy /mZ(co) = 
aco 2 with a = 0,0.2,0.4,0.6,0.8 and 1.0, respectively (the K-K re- 
lated ReY,(a>) is included). A constant DOS ATq(co) = 1.0 was chosen 
for the non-interacting case (a = 0). (b) The results of calculations 
of — i Y,i c Im{G(k, co)(l — with the same self-energies as in (a). 



contribution of all closed-linked skeleton diagrams, 
but with replacing all Green's function lines by 
the full Green's functions G r (co„). 

(8) 

The explicit expression of £2 was given in LW(50) (refers 
Eq.(50) of Ref.[2]). The functional £2 was ingeniously de- 
signed to satisfy the famous variational theorem of the free 
energy functional: 



3D. 



0. 



(9) 



And this theorem can be easily proven if the functional £2 
satisfies the following variational property 



do! 



3E r (co„) 



r£[G,-(co„)] 2 E r (co„) 



(10) 
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which was shown in LW(51). Using the same reasoning used 
to derive Eq.dTOb. we can derive a similar variational relation 



because the leading temperature dependent parts in the re- 



3£2 
3;C0„ 



-r^[G,-(co„)] 2 (l-^-)E r (co„). (11) 



maining terms —TY, rn e"°" [G r (co„)E,-(co„)] + £2 cancels ex- 
actly each other. Hence, S(T) considered by Luttinger is the 
following 



The above relation will be used for our purpose later and the SLuttinger(T) — 

proof of it is trivial if we note the expression of G,7 1 = /co„ — 

e, - E,.(co„) and replace ^- in Eq.(lO) by = 

In order to calculate the entropy from the free energy func- 
tional Eq.(|7]l, we need to extract the leading temperature de- 
pendent parts of it. Using a standard method of the Matsubara 
frequency summation, Eq.(Q is written as 
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-G r (z)Mz)}-f^n F (z)a'(z) (12) 



where all Matsubara frequencies are analytically transformed 
to complex numbers as /co„ = z and the functional £2 (z) is 
also understood as £2 (/co„ — > z) replacing the overall Mat- 
subara frequency summation —T^ n in the original definition 
of the functional £2 (/co„) by the contour integral § |I;«f(z) 
with the Fermi-Dirac distribution function n F (z). Now it is 
clear that there are only two places where contain the temper- 
ature dependence in the above free energy functional Eq.(IT2b: 
rif(z) and E r (z). As Luttinger argued [3], the leading temper- 
ature dependence should come from the explicit summation 
of /CO,, (equivalently in n F (z)) and the temperature variation 
of "E,-(r) — E, .(T = 0)" is a higher order and should be ne- 
glected. Therefore, using S(T) = —dQ.(T)/dT, we can write 
down S(T) as follows 



S(T) 
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£/w{ln G7 1 (co) + G r (co)E° (co) } 
7m£2' (co) , 
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(13) 



where the contour path of § is deformed along the real fre- 
quency axis a la the appendix of Ref.[2]; the CO integration for 
[— oo,oo] should be carried infmitesimally above the real axis, 
i.e. for co + /r|. E°(co) means E r (co, T = 0) and it is under- 
stood that every E, (co,r) implicit in the above expression is 
replaced by E'?(co). While the above equation for S(T) is the 
undoubtedly correct expression, Luttinger argued in Ref.[3] 
that the leading temperature dependence of £2(T) (Eq.©) is 
contained only in 



£2 



temp 



(T) « -r^ e im " 0+ ln[e r + E r (co„) - /co„] (14) 



Note the difference of the expressions of the entropy S(T) 
(Eq.dOll) derived in this paper and the one SLuttinger(T) 
(Eq.dT5t) considered by Luttinger in Ref.[3]. Below we will 
show that the exact cancellation of the remaining terms - 
claimed by Luttinger - is incomplete and therefore we have 
to keep the full expression of the entropy Eq.([T3l. Expectedly 
the calculation results of the SH coefficient y from Eq.([T3l 
and yiuttinger from Eq.([T5]) are totally different: the former one 
gives y unrenormalized regardless of the strength of the inter- 
action while the latter one gives an enhanced ytuttinger pro- 
portional to the interaction strength as widely believed in the 
community ever since the proof of Luttinger 13|] . 

To obtain y= lim r ^oC(r)/r = Hm T ^odS(T)/dT, we 
only need to extract T-linear terms in S(T) or S Luttinger {T). 
Utilizing Sommerfeld expansion, we then only need to ex- 
tract co-linear terms in the integrand of Im... in S(T) or 

S Luttinger (T). Let US first calculate y Luttinger from S Luttinger- The 

leading Taylor expansion of the integrand in SLuttinger is 



7m{ln[e r +£?-<»]} = 7mG r (co = 0)(l 
+ O(co 2 )... 



a/^E° 



3co 



(0=0 



)co 
(16) 



Using the notation ( 1 



C0=0 



) = 1 + X r = Z r , we obtained 



--TlmG r (<o = 0)(l+X r )ttZ-N(0)=Z-No(0) (17) 

71 r 

where Z is a Fermi surface average of Z,-. In Fig. 1(b), we 
showed the numerical calculations of — ^^ r /m{G r (co)(l — 

-^)} with the varying interaction strength. This quantity has 
no direct physical meaning (it becomes even negative at higher 
energies) but its zero frequency value NLuttinger{(£> — 0) clearly 
demonstrated the result of Eq.(17) and showed what quantity 
was used by Luttinger for the calculation of the SH coefficient. 
Substituting the result of Eq.(17) into Eq.(15), we can derive 
the same result as Luttinger had obtainedyl] as 



y Luttinger = —Z-No(0), 



(18) 



so that the SH coefficient yLuttinger is indeed enhanced by the 
factor Z compared to the non-interacting case. 

Now let us use the exact expression S(T) of Eq.(13) to de- 
rive y. The coefficients of the co-linear terms of the integrand 
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of Eq.(13), {lnG^co) + G r (co)E?(co)} -£l'(co), are the fol- 
lowing 

3E° 



[GJ 2 «d-^) 



3co 



(19) 



Above we have arranged the expansions of each three terms 
lnG^Co), G r ((o)E°(co) and — Q. (co) into three separate lines 
for clarity. There are lots of cancellations and the final re- 
sult should be compared to Eq.(16) obtained from Si utt i nser (T) 
Eq.(15). In fact, the above cancellation is the consistent result 
of the variational theorem of Eq.(9) which requires that all 
variations of 3E r in the total free energy functional £2 should 
sum up to zero yfl. In this sense, the expression of SiMtinger 
in Eq.(15) with Eq.(16) cannot be correct since it contains 3E r 
term. 

Now it is a trivial matter to calculate the SH coefficient y 
substituting the result of Eq.(19) into Eq.(13) as 



-7mG,-(co = 0) 
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-■N(0) = — -iVo(O) 



(20) 



where we have counted the spin summation factor 2 in the 
index r = (k,a). The above result shows that the SH coeffi- 
cient y of the interacting Fermi system is proportional to the 
interacting DOS at the chemical potential N(0), which is con- 
sistent with our physical intuition. However, due to the ab- 
sence of renormalization of the DOS at chemical potential, 
i.e., N(0) = No{0), which we have shown in Eq.(5), the SH 
coefficient y is not renormalized by the interaction as far as 
the system remains a Fermi liquid state. 

Other Physical Quantities - The renormalized q.p. mass 
m* pa Zmo due to interaction is measured by different experi- 
mental probes. Indeed the energy dispersion of the q.p. pole 
E(k), defined by CD — e(k) — E(&,C0) = 0, is renormalized as 
E(k) rj z{k)/Z and should be directly measured by ARPES 
without any interpretation or confusion. Another common 
tool to measure m* is the dHvA effect with the applied exter- 
nal field H. In this case, the effective mass is measured from 
the temperature reduction factor of the signal strength which 
is given by the Lifshitz-Kosevich formula Rj- ~ exp (— T / co c ) 
[Q], where C0 C is the cyclotron frequency. co t is determined 
by the q.p. energy distance between the Landau levels quan- 
tized by the field H as AE = h03 c , and the Landau level is 



determined by the q.p. dispersion E(k) to the first approx- 
imation, hence co c = eH/m*c. Therefore, the dHvA effect 
measurement can provide an information of m*. Lastly, the 
optical spectroscopy measurements need a more careful inter- 
pretation. The total spectral density near the Fermi level is 
not enhanced by interaction as shown in Eq.(5), however the 
width of the q.p. dispersion is narrowed by the factor 1 /Z as 
shown in Fig. 1(a). Therefore, for example, the width of the 
Drude spectra in the optical conductivity should be reduced 
by the factor 1/Z , while the absolute magnitude of the zero 
frequency conductivity o(co = 0) nor the total Drude spectral 
weight is not expected to be enhanced. However, because the 
optical conductivity is a transport property, it is necessary also 
to count on the scattering rate due to the interaction besides the 
q.p. DOS. Therefore, for more complete details of the optical 
properties of the interacting fermion systems, we need to an- 
alyze the two particle correlation function which was beyond 
the scope of the current paper. 



Conclusions - In this paper, we have proven that the SH 
coefficient of the interacting fermion system is not renormal- 
ized by the interaction and remains basically the same value 
as the non-interacting limit of the system. This result is in 
stark contrast to the longtime accepted idea of the interaction- 
enhanced SH coefficient since the proof of Luttinger in 1960 
J^t]. The implications of our finding should be far reaching 
and profound because the enhanced SH coefficient y in the in- 
teracting Fermion systems has been accepted and utilized for 
the last 50 years as a pivotal building concept in the study of 
the interacting Fermi liquid systems both in theory and in ex- 
periment. We need to rethink many of the previous ideas and 
measurements based on this -now proven wrong - concept. 
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